Scaling for domain growth in the Ising model with competing dynamics 
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Ch '. 

We study the domain growth of the one-dimensional kinetic Ising model under the competing 
influence of Glauber dynamics at temperature T and Kawasaki dynamics with a configuration- 
QO ' independent rate. The scaling of the structure factor is shown to have the form for nonconserved 

dynamics with the corrections arising from the spin-exchange process, i.e., S(k,t) — Lgo(kL,t/r) + 
gi(kL,t/r) + and the corresponding scaling functions are calculated analytically. A correction 
to the Porod law at zero temperature is also given. 
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The kinetic Ising models have been widely studied to understand the far-from-equilibrium phenomena such as the 
nonconserved Glauber model jlj] and the conserved Kawasaki model Q. In recent years, the systems with more complex 
mechanisms, e.g., the Ising models with competing Glauber and Kawasaki dynamics HQ, have been introduced to 
investigate the basic questions of nonequilibrium phase transitions and critical phenomena. For these models with 
competing dynamics, much of the interest has been focused on the properties of nonequilibrium steady state, including 
the nonequilibrium phase diagrams and the critical exponents HQ- However, the problem of domain growth, 
which describes the behavior of the system quenched into an ordered phase from a high-temperature initial state and 
is important in understanding the dynamics of nonequilibrium processes, has received relatively little attention. 

Our interest here is in the study of domain growth for the Ising system with competing dynamics. It is well known 
HI that in the late stage regime, the nonequilibrium process of domain growth exhibits dynamic scaling behaviors 
and the scaling forms of the equal-time pair correlation function C (r, t) and its Fourier transformation, the structure 
' factor S(k,t), are given by C(r,t) = f{r/L{t)) and S(k,t) = L(t) d g (kL(t)), respectively, where L{t) is the single 
, length scale characterizing the domain structure and d is the spatial dimensionality. The corrections to the above 
scaling arising from the nonscaling initial condition were determined recently at zero temperature j[o|. However, 
the direct demonstration of scaling and the exact calculation of the scaling functions are still lacking except in some 
simple models fil]|-|l3| . Thus it is of interest to give a direct and exact calculation of scaling behaviors for the system 
with more complex dynamics. 

The system considered here is the one-dimensional Ising model evolving by a combination of the spin-flip Glauber 
process at temperature T and the spin-exchange Kawasaki process with a configuration-independent rate that occurs 
as if the system were at infinite temperature ||,|l4| . In this work we study the domain growth of this soluble model 
following a sudden quench from an initial high temperature to a final temperature T and calculate the corresponding 
scaling properties analytically. 

The Hamiltonian for this one-dimensional system is given by 
' 

Ch ' N 

O'. H = -jJ2<n°i+i, (!) 



where <7j = ±1 is the Ising spin variable. The configuration of the system evolves with time via the combination of 
Glauber dynamics with the spin-flip rate = (2ti)~ 1 [1 — 70^ (<Tj_i + Ci+i) /2], where 7 = tanh2ii~ with K — 

J/krjT, and Kawasaki dynamics with the configuration-independent spin-exchange rate W^li = (2t2)~ 1 (1 — OjCTj+i) . 
The equations for the expectation value of the spin and the equal-time pair correlation function have been derived 
fhij and the latter is written as 

T'dC{r, t)/dt = -2C(r, t) + 7' [C(r - 1, t) + C(r + 1, *)] (r > 2), 
r'dC(l,t)/dt = ~^C{1, t) + j'C(2, *) + y - 

1 + K 1 + K 

C(0,i) = l (2) 
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provided the initial probability distribution is translationally invariant, where the pair correlation function C(r, t) = 
(a t (t)a l+r (t)), K = 2ri/r 2 , r' = n/(l + k), and f = (7 + k)/(1 + k). 
From Eq. (Q) the scaling results for domain growth can be calculated directly following the method given by Bray 
By Fourier transforming in space, the equation for the structure factor S(k,t) is obtained, 



r'dSik, t) Idt = -7fc5(fc, t) + A(t) n(t) cos k 

1 + K 



(3) 



where 7^. = 2(1 — 7'cosfc), A(t) = ^2i c 'ykS(k,t)/N, and n(t) = [1 — C(l,t)]/2 denotes the average wall density. 
This equation can be solved by a Laplace transformation in time via S(k,s) — dtexp(—st)S(k,t). In the scaling 
regime of late times after the quench, the characteristic domain size L(t) is known to obey the growth law L(t) ~ t x , 
where the growth exponent x is 1/2 and 1/3 for systems with nonconserved and conserved scalar order parameters, 
respectively H[§. Its inverse, the average wall density, then satisfies n(t) = bt~ x /2 = bt^^ 1 /2 with < v < 1, as well 
as a certain coefficient b, and the corresponding Laplace transform is n(s) = bT{v)s~ l/ /2. Here we assume that the 
above scaling forms of the domain size and the wall density remain for this model with competing dynamics, which 
will be justified a posteriori. Thus, on the condition that the initial state contains no long-range order, in the scaling 
limit (s — > 0, k — > with sjk x l x arbitrary) we have 
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when domains grow at small but nonzero temperature T. Here k! = ~ 2 exp(— 2K) / y/1 + k, where £ is the 
correlation length for this system with competing dynamics fl4|| . 
Using the inverse Laplace transform on Eq. (H) yields 
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where B(x,y) is the Beta function and \Fi{a\ /?; z) is the degenerate hypergeometric function Jlq ). It is shown from 
Eq. (||) that a scaling solution requires L(t) ~ i 1 / 2 , i.e., = 1/2, and then the equal-time structure factor has the 
scaling form 



S{k, t) = t x l 2 g a {k 2 t, K l2 t) + gi{k 2 t, k'H) + r 1 l 2 g 2 (k 2 t, n' 2 t) + t^g^kH, k'H), 
with the scaling functions 
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g 3 (k 2 t, k'H) 
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Thus, from this scaling result of the structure factor it is shown that the domain size L(t) scales as t x and, corre- 
spondingly, the wall density n(t) has the scaling form t~ x , which is consistent with the assumption made above before 
the calculation. Here the growth exponent x is 1/2 for the competing dynamics, the same as that for the (single) 
nonconserved Glauber system. 

For K l2 t 3> 1, one can obtain the asymptotic result of S(k,t), that is, 
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while for K l2 t — > 0, the scaling result for quench to T = 0, which is the main interest in the one-dimensional system, 
is given by 
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which is of the scaling form 
For k 2 t > 1, Eq. (|l|) becomes 
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which yields a correction to the Porod law |16| . 

It should be noted that the scaling forms of the structure factor (|^) at nonzero temperature and Eq. (|l3|) at zero 
temperature can be rewritten as 

S(k,t) = Lg (kL,t/T) + 9l {kL,t/T) + L- 1 g 2 {kL,t/T) + L~ 2 g 3 (kL,t/T), (15) 

with r = 1/k' 2 = £ 2 and L - t 1/2 , and 

S(k, t) = Lg Q (kL) + 9l (kL) + L- 1 g 2 {kL), (16) 

respectively. By Fourier transforming, the corresponding scaling of the equal-time pair correlation function is obtained 

C(r, t) = f (r/L, t/r) + L^f^r/L, t/r) + L- 2 f 2 (r/L, t/r) + L~ 3 f 3 (r/L, t/r), (17) 

at nonzero temperature and 

C(r,t) = f (r/L) + L^Mr/L) + L- 2 f 2 (r/L) (18) 

at T = 0. In the expressions of these scaling forms, the first leading term exhibits the scalings for the nonconserved 
dynamics in one dimension |9|Jl2]| associated with the spin-flip Glauber process, while the other terms are the correc- 
tions to scaling that arise from the spin-exchange Kawasaki process, which can be shown from the expressions of the 
scaling functions. In formulas (|§|)-(|lO|) for quench to nonzero temperature, the correction functions gi, g 2 , and g 3 vary 
directly as k/(1 + k) and vanish if the contribution of the spin-exchange process is negligible or, equivalently, k — > 0. 
Consequently, the scaling form S(k, t) — t 1 ^ 2 go(k 2 t, K ,2 t) with the scaling function go(x, y) and the asymptotic result 
S(k,t) — > 2k' /(k 2 + k' 2 ) for K ,2 t ^s> 1 obtained by Bray |Q for the one-dimensional Ising model with the (single) 
Glauber dynamics are recovered from Eqs. (|6|)-(pd[). Similar results can be obtained for zero temperature, that is, one 
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can recover the well known scaling form S(k, t) — t 1 / 2 'qp(k 2 t) and the Porod law S(k,t) ~ k~ 2 for Glauber Ising chain 
|^2| , ^3| from Eqs. (p^|)-(14). Moreover, from Eqs. (|l7| ) and ( |l8| ) we note that the leading correction to the noncon- 
served scaling form of the pair correlation function is L _1 /i(r/L) here, originating from the spin-exchange process, 
while in the absence of this process, associated with departure of the initial condition from the scaling morphology, it 
is L~ 4 /i(r/L), according to the recent work of Bray et al. Go ]. 

From the above scaling solution no correction to scaling of the domain size L(t) and the wall density n(t) is shown, 
as expected before the calculation. However, if assuming the correction to scaling of L(t) and n(t) before solving the 
evolution equation (||), e.g., n(t) ~ b\t~ Xl + b^t^** + ■ ■ ■ with < x\ < X2 < ■ ■ •, we find that the scaling solution 
similar to Eqs. (^J) and ( |l5| ) cannot be obtained. 

In summary, we have shown that for this Ising system with competing Glauber and Kawasaki dynamics, the spin- 
exchange process with a configuration-independent rate has the effect of corrections to scaling, while the leading 
contribution to scaling is of the form for the nonconserved Glauber dynamics, as shown in Eqs. (|^) and (|l2|). More 
work is needed for a further understanding of these results. 
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